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Abstract. We illustrate how to compute asymptotic interactions between 
discrete solitary waves of dispersive equations, using the approach proposed by 
Manton [Nucl. Phys. B 150, 397 (1979)]. We also discuss the complications 
arising due to discreteness and showcase the application of the method in nonlinear 
Schrodinger, as well as in Klein-Gordon lattices, finding excellent agreement with 
direct numerical computations. 
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1. Introduction 

Solitary wave interactions and collisions are among the signature characteristics of such 
nonlinear coherent structures and have been extensively studied in both integrable 
and non-integrable dispersive nonlinear models [fl[21 ED ED HI EH • The relevant physical 
contexts where such phenomena may arise range from the more traditional areas of 
fluid mechanics and plasma physics (discussed in the above references) or high-energy 
physics [7| to nonlinear optics [SJ , and the recently emerging applications of Bose- 
Einstein condensates [Q3 [TJ . 

On the other hand, in the past decade there has been an explosive parallel 
evolution of the theory and applications of discrete nonlinear waves, namely discrete 
solitons and discrete breathers, see, e.g. |12j for a number of review papers on the 
topic. In this case, the applications refer to nonlinear waveguide arrays [Qj and 
photorefractive crystals [Q] , to Bose-Einstein condensates in deep optical lattices [T^] , 
to the local denaturation of the DNA double strand [T^j, and to micro- mechanical 
cantilever arrays [T7] among many others. 

Our aim in the present work is to illustrate a technique for computing solitary 
wave interactions in discrete systems. We should mention here that for continuum 
systems, there exists a large variety of methods for computing such interactions. 
These range from the perturbation theoretical works of , to the variational methods 
of [191 121)1 |2~T] . the Fredholm-alternative based technique of [221 or the more rigorous 
calculations of [23] based on Lin's method. In a recent publication |24|. we took 
an alternative route to these methods, by implementing an asymptotic calculation 
using the approach proposed by Manton |25j (which, in turn, was generalizing the 
earlier work of |2fi|). Performing such calculations in discrete systems is, however, 
considerably less straightforward. One of the complications is that, aside from the 
tail-tail interaction of the waves, each wave may also reside on a periodic substrate 
potential, usually termed the Peierls-Nabarro barrier [T5|. As a result, except for the 
(typically exponential in the wave separation) tail-induced potential, a local potential 
may arise due to discreteness, resulting in a washboard potential structure which can 
be captured by variational techniques [27)1 l2"T] , 

In what follows, we show how to compute discrete solitary wave interactions by 
means of Manton's method. We bypass the issue of the Peierls-Nabarro barrier by 
working with problems whose static solutions do not encounter such a barrier (so 
that they can be placed anywhere on the lattice). We demonstrate this approach 
both for breathers in nonlinear Schrodinger lattices - namely in its famous Ablowitz- 
Ladik (ALNLS) variant [22] - and for kinks in Klein-Gordon lattices - namely in a 
discrete version of the <fi 4 model, proposed in [2B1 [213 EH- While the former model 
is integrable, integrability is not a key ingredient in our calculation. Instead, the 
existence of a discrete analog of momentum conservation law (and also the absence 
of the Peierls-Nabarro barrier) are indispensable for being able to carry through the 
Manton calculation as we show in the following sections. We compare our analytical 
results with direct numerical computations, finding excellent agreement between the 
two during the time interval in which the waves interact without losing their individual 
character. Furthermore, our calculations also capture the correct continuum limit 
(known from the earlier works of |19l I25j : see also [211), as the lattice spacing h — > 0. 

Our presentation is organized as follows. In section 2, we illustrate the method 
for the breathing solitons of the Ablowitz-Ladik model. In section 3, we compute 
the interaction potential for kinks in a non-intcgrablc, discrete <fi A model. Finally, in 
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section 4 we summarize our results and present our conclusions and some interesting 
directions for future study. 

2. Soliton Interactions in the Ablowitz-Ladik Model 

Similarly to the corresponding continuum calculation (see, e.g., |24| for such 
calculations in different models), the key to performing a calculation for the discrete 
case using Manton's method is the existence of a momentum as a conserved quantity. 
In the ALNLS model, such a momentum operator is one of the integrals of motion 
(that is why integrability is desirable for performing such a calculation, even though, as 
we will show below, it is not a necessary condition), and has the form (see, e.g., |2"ll31| 
for details) 

oo 

P = i hPn^n+l ~ Ipn^n+l) , (1) 

n— — oo 

where ip is the complex field, and the star is used to denote complex conjugation. We 
will focus here on the ALNLS equation in the form [3*T] 

#n + (V , n+l+Vn-l)(l+IV'n| 2 ) =0, (2) 

with the explicit stationary soliton solution 

ijj n = sinh(/?)sech ((3(n — s j) exp(ier), (3) 

where a — 2 cosh(/?). Here the overdot denotes time derivative, (3 is the inverse width 
of the pulse soliton and s is the arbitrary position of the center of the pulse. We 
commence by examining in-phase solitons, but we will also relax this constraint later. 
Note that Eq. J5J has the lattice spacing h scaled out, but it is straightforward to 
incorporate it and we will see how to relate it to the relevant continuum limit of 
h — > for our calculations. 

In the spirit of [35], we now consider two solitons, one centered at and 
one centered at s 0, i.e., two widely-separated solitons. We compute dP/dt by 
performing the summation over n not for the infinite lattice (when the result would 
be zero due to the relevant conservation law), but rather from n = M to n = N, with 
M < 0, and ^ N s. The idea behind this calculation is that, in fact, the force 
in this interval is not going to be zero, but rather would be finite due to the soliton- 
soliton interaction. When summing on the infinite line, as a result of Newton's third 
law, the action of the first soliton on the second and the equal and opposite reaction 
of the second on the first cancel each other, thus resulting in zero net momentum gain. 
However, for a finite interval encompassing only one soliton, the amount of momentum 
gain is finite, due to the fact that the one soliton experiences the pull (or push) of 
the other soliton at the boundary of the interval where we perform the calculation. In 
precise mathematical terms, we evaluate: 

rW N 

^= -2£ (l^ + i| 2 -hM 2 ) 

n=M 

N 

+ + V^«+2) (1 + IV-n+ll 2 ) 

n=M 
N 

- (^n-lK + l + r n -li>n+l) (1 + l^nl 2 ) • (4) 

7l=M 
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However, observing the telescopic nature of the sums in the right hand side (RHS) of 
Eq. (0J, we infer that 

= -2(|^ +1 | 2 -hM 2 ) 

+ (^N^N+2 + 1pN^N+2) (l + |V^V+l| 2 ) 

- (l/'M-lV'M+l + 1>M-l1>M+l) (1 + I^A/I 2 ) ■ (5) 

As usual in Manton's method, and based on intuitive physical arguments, the main 
contribution in this asymptotic calculation stems from the boundary between the 
two solitons. Hence, we drop the terms with subscript M and only consider the 
contributions with subscript N in what follows. 
We then invoke the soliton ansatz 

i> n = ^ + ^ (6) 

with ipn^ = sinh(/?)sech(/3n) exp(za) and — sinh(/?)sech(/3(n — s)) exp(ier) (i.e., 
two in-phase solitons). Since <C N <C s, we can use the asymptotic form of the 
soliton tail at n = N, according to 

ipW = 2 sinh(/3) exp(-/3AT) exp(ier) , (7) 

</>i 2) = 2 sinh(/3) exp(/3(iV - s)) exp(icr) . (8) 

Substituting the ansatz of Eq. © and the expressions in Eqs. 10-© into Eq. JSJ), 
and computing the terms arising from the soliton-soliton interaction, we obtain that 
HP 

— « 32 sinh 4 (/3) exp(-/3s) . (9) 

We note a number of comments on this calculation. 

• Equation of motion for s(t): As explained in |24j . in order to obtain the equation 
of motion for the inter-soliton separation s(t), we can use Newton's equation in 
the form 

dP , . 

M s s = -2-, (10) 

where M s is the mass of the soliton; the factor "2" comes from the fact that 
there is an equal and opposite pull (or push) on the second soliton, and hence 
their relative distance decreases by twice the contribution of dP/dt to each of 
them; and finally the unfamiliar "-" sign originates from the fact that a positive 
boundary contribution to dP j dt decreases the soliton distance, while the opposite 
is true for a negative dP j dt. In the present case, the soliton mass is given by 

oo 

M s = ln(l + |V^| 2 ) -2/3 (11) 

n— — oo 

(once again, a direct benefit of integrability being the immediate accessibility of 
the relevant conservation law). Then one can straightforwardly infer the equation 
for s(t) as 

s = --sinh 4 (/?)exp(-/? S ), (12) 

while the relevant effective soliton interaction potential (for a unit mass particle) 
is 



V(s) = --2 S mh i ([3)eM-Ps)- (13) 

P 
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Figure 1. The top panel shows the soliton-soliton separation as a function 
of time, s(t), for the numerical integration of the ALNLS equation with two 
superposed solitons (at n = —5 and n = 5) as the initial condition. The solid line 
yields the numerical result (calculated by obtaining independently the center of 
mass of the left and right pulse in the configuration, and subtracting the former 
from the latter), while the dashed line provides our asymptotic prediction of 
Eq. 1121 . Note the excellent agreement between the two, essentially up to the 
collision time, since the two curves cannot be distinguished, practically, up to 
that time. The bottom panel shows the space-time contour plot of the "local 
mass" m n = ln(l + \ip„\ 2 ), clearly illustrating the evolution towards the collision 
of the center of mass of the two solitons. 



We have examined the validity of the pertinent calculation by studying 
numerically the two-soliton collisions in the ALNLS model. Typical results are 
illustrated in Fig. \T\ior (3 = 2 and two solitons initialized at a distance of s(0) = 10 
between them. The upper panel of the figure shows the numerical (solid line) 
versus theoretical (dashed line) prediction for the soliton separation as a function 
of time, s(t), revealing excellent agreement for all times until about t « 657, when 
the separation becomes s < 3 and the solitons can no longer be characterized as 
individual entities; hence, it is reasonable that our asymptotic calculations fail at 
that point. In fact, it appears that the analytical calculation gives a very good 
agreement with the numerics well past the point where we might have anticipated 
such agreement (on the basis of our expansion assumptions), and approximately 
even up to the collision point. 

• Role of the soliton relative phase: In order to examine it, it is simple to impart 
to one of the ALNLS solitons a free phase exp(i</>), with respect to the other 
soliton. In this case, the results of Eqs. H12(l - I|13(l are retrieved, but with the 
RHS expression multiplied by cos(0). This is rather natural as in-phase ALNLS 
solitons are expected to attract each other, while out-of-phase ones are expected 
to experience mutual repulsion |191 1201 |2~T] . 

• Retrieving the continuum limit: In the continuum limit, one should recover the 
relevant expression of |32| (see also |24|L In reshaping our expression to match 
the latter, we observe that in the continuum limit the discrete spacing s will 
appropriately renormalize to the continuum spacing, while sinh(/3) — * (3 for 
the continuum soliton. Finally, our prefactor of 32 in Eq. will match the 
corresponding prefactor of 8 in |2 II I32| , since the momentum is defined with a 
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factor of 1/2 (as in Eq. (11) of [21]) and there is an extra factor of 1/2 multiplying 
the RHS of Eq. @ to match it with the (continuum limit of the) RHS of the 
equations used in |2"31I32"] . 



3. Kink Interactions in a Discrete <j> Model 

We consider now a rather different class of models, namely Klein-Gordon lattices. In 
particular, we focus on a discretization of the <p A field theory of the form: 

il n = CA 2 u n + 2u n - u 2 n (u n+ i + u n -x), (14) 

where C = 1/h 2 is the coupling parameter and A2U„ = (u n+ i + it n _i — 1u n ) stands 
for the discrete Laplacian. This type of discretization for the RHS was discussed 
in [2BI (motivated by ODE examples), and independently rediscovered for the </> 4 field 
theory in |29II3U| . An interesting feature of the model is that, while non-integrable, it 
possesses exact solitary wave solutions in the form 

u n = tanh(a(n — s)), (15) 

where s denotes the center of the kink, a is the inverse width of the kink and 
cosh(2a) = (1 + h 2 ) /(l — h 2 ). This feature, while useful, is again not an indispensable 
one for the development of the Manton procedure. In fact, the only related property 
that is necessary for the calculation is the exponential tail of the solitary waves (which 
can, in fact, be derived by an appropriate exponential tail ansatz, even if it is not 
explicitly available in the form of an analytical solution). 

Another property of this model, more crucial for our considerations, is that 
it possesses a discrete analog of the continuum momentum conservation law. In 
particular, as shown in |29j . the relevant momentum is of the form: 

^ oc 

P= ~n X! " n _ U n-l) ■ ( 16 ) 

n— — oo 

One can then perform a similar calculation of dP/dt summing from n = M to n = N, 
thus obtaining: 

dP 1 N 

= >^ [H(u n+ i, u n ) - H(u n , w„-i) + u n u n+ i - u n -iu n ] (17) 

n=M 

where 

H(u n+1 ,u n ) = C {u n+ i - u n ) 2 + 2u n u n+1 - u 2 n u 2 n+l . (18) 
Given the telescopic nature of the summation, we can easily infer that 
dP 1 

— = -- [H(un+i,u n ) - H(um,u M -i) + unu n +i - um-\Um] -(19) 
We now use the kink-antikink ansatz of the form (see also |24l 125) ) 

u n =uW + u$ +r, (20) 

where r is the value of the inhomogeneous background steady state on which the kinks 
exist. In the case of the cf) 4 model, it can be easily inferred that r = ±1. Without loss 
of generality we use r = — 1 here. The asymptotic form of the kink profiles, using the 
expression of Eq. (|15|l at n = N, will be: 

(1) - 1 - 2exp(-2aA0 , (21) 



'N 

{ (2) 



1 -2exp(2a(AT-s)) . (22) 
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Notice that in general the first term in the RHS of Eqs. (|21|) - (|22[1 is — r. 

Using the ansatz of Eq. i|2(J|) along with Eqs. I|21|) - (|22ll in the expression of Eq. 
11911 , and isolating the contributions stemming from the tail-tail interaction (to leading 
order, once again coming solely from n — N), we obtain that 

dP 

— « 16cxp(-2as) [Csinh 2 (a) + cosh 2 (a)] . (23) 
We again make a number of comments on the calculation. 

• Equation of motion for s{t): In this model, we do not have an immediately 
available definition of the soliton mass M s . However, considering the momentum 
itself, we realize that using a time dependent center of mass s(t) for the position 
of the center of one of the solitons, and substituting that along with Eq. I|15|) in 
Eq. (|16J) . we obtain 

oo 

P = as sech 2 (a(n — s)) [tanh (a(n + 1 — s)) — tanh (a(n — s))] . (24) 

n— — oo 

Since a is an effective spacing parameter that will be used to renormalize the 
distance in the continuum limit (see below), we will consider the expression 
multiplying as in the above formula of the momentum as the mass of the kink 
Mfc. Mk and a have been computed numerically for various spacings h in the 
interval (0, 1) and are shown in Fig. [21 It can be clearly observed that Mk tends 
to its continuum limit of 4/3 (= u^dx = scch(x) 4 da;), validating the 
usefulness of this definition. Note also in the same graph that the behavior of 
a is linear in h for small enough h (in fact, the relevant Taylor expansion is 
a w h + h 3 /3 + . . . for small h). Let us also mention that while the expression of 
Eq. (|15fl for the exact solution appears to be necessary for the computation of the 
kink mass from Eq. I|24|) . this constraint can be relaxed as well. For instance, one 
can quite efficiently approximate the mass numerically (see also our numerical 
calculations below). 

Equating the right hand side expression from the momentum of Eq. I)23|) [doubled 
and with a minus sign for reasons similar to those mentioned previously for Eq. 
(|l(J|l ]. with the expression stemming from the time derivative of P = asMk, we 
obtain the evolution dynamics 

32exp(-2as) r ^ . , 2/ % , 2, m ,^^\ 

s = — rCsmli (a) + cosh 2 (a) . (25) 

M k a 1 1 

The relevant effective tail-tail kink-antikink interaction potential (for a unit mass 
particle) can then be obtained in the form 

= _ 16exp(-2a£) r Csinh 2 (a) + cosh 2 (a) n _ (2g) 

These predictions are tested against numerical simulations, obtaining once again 
excellent agreement almost up to the collision point (the two trajectories only 
separate visibly for t w 45, when the distance of the kink and anti-kink is less 
than 5 sites). The results are illustrated in Fig. [21 

• Retrieving the continuum limit: In the continuum limit, we have to compare 
Eq. ([23 with the expression for dP/dt that was obtained in {25) . 

dP 

— w 2A 2 m 2 exp(-ms) , (27) 
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Figure 2. The graph shows the dependence of the kink mass given by Eq. 
1241 (solid line) as a function of the lattice spacing h. The continuum limit 
corresponding to Mj. =4/3 (see the main text) is shown by the dash-dotted line. 
The dashed line illustrates the dependence of the parameter a on the spacing h, 
while the dotted line shows the linear behavior for comparison. 
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Figure 3. Same as Figure 1, but now for the kink-antikink interaction of the </< 4 
model for h = 0.5. The bottom panel shows the space-time contour plots of the 
quantity D n = Ji n +i — u n -\ whose square we used for our local (approximate) 
calculation of the centers of mass of the two kink and antikink, respectively. 






for two solitons with the asymptotic form 

= - r + Aex-p(-mx), (28) 

u (2) = - r + A cx-p(m(x - s)). (29) 

Hence, when renormalizing the distance x by our a close to the continuum limit 
of h — > 0, and with m — 2 and A — —2, one obtains dP/dt w 32exp(— 2s) 
for the <p A field theory. We retrieve this result from the expression of Eq. 
since lim/^o C(sinh(a)) 2 = 1 and cosh(2a) — > 1 as h — * 0, by a straighforward 
limit in the expression defining a (just below Eq. I|15|l 'l. Once again, the correct 
continuum limit is obtained from the discrete expression as a special case. 
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4. Conclusions 

In this work, we have presented a self-contained description of the way to use the 
technique first proposed by Manton in [23 for the continuum Klein-Gordon models, 
in order to quantify tail-tail interactions of soliton solutions of discrete, dispersive 
nonlinear wave equations. We demonstrated the calculation via two prototypical 
examples, one in the form of soliton-soliton interactions in the completely integrable 
Ablowitz-Ladik model and one in the form of kink-antikink interactions in a (non- 
integrable) Klein-Gordon chain. We highlighted the key necessary ingredients for 
carrying out such a calculation, namely the existence of a momentum-like quantity 
which is conserved due to its "telescopic property" . This allows, similarly to the 
continuum limit, to arrive at an expression involving only boundary terms when 
computing the force exerted on a fraction of the lattice encompassing one of the two 
solitary waves. Assessing the leading order contribution from these surface terms and 
equating it (with appropriate prefactors taking into consideration the total magnitude 
and direction of the relative change of displacement of the two waves) to the particle- 
like momentum of the coherent structures, we obtain the kinematics of their relative 
displacement. The relevant expression clearly highlights the exponential nature of the 
tail-tail interactions (which is natural, given the exponential tails of the solitary waves). 
The obtained expressions have been tested both against the limiting case of their 
continuum counterparts which were previously available, and against direct numerical 
simulations, providing excellent agreement with the direct numerical experiments, 
practically up to the wave collision point (for the attractive interactions considered 
herein). 

The success of our predictions against the corresponding numerical simulations 
naturally raises the question of generalizations, as well as limitations of the technique. 
The most significant limitation is that "standard" discretizations often defy the 
existence of a momentum conservation law, and do not share such a property with their 
continuum siblings, restraining themselves to the integer shift invariance, rather than 
an effective continuum-like translational invariance. It then takes either integrable 
(such as e.g. the ones of [22]) or "non-standard" (see e.g. [2H1 CS3 and references 
therein) discretizations to circumvent this difficulty and possess a discrete analog of the 
momentum. Hence, it would be very desirable to be able to extend our considerations 
to cases where the momentum conservation is absent, hopefully obtaining, in addition 
to the tail-tail terms already captured above, the ones from the local, periodic Peierls- 
Nabarro barrier. Another natural extension of the Manton approach (in fact, both 
in the continuum and in the discrete systems) would be the study of its multi- 
dimensional analog. While by no means straightforward (especially given the sparsity 
of analytically available solutions and the difficulty of questions such as the definition of 
the right momentum/contour), it constitutes an increasingly relevant generalization of 
the considerations presented herein. Work along these directions is currently underway 
and will be reported in future publications. 
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